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Abstract 

C 4 photosynthesis evolved independently numerous times, probably in response to declining atmospheric C0 2 con- 
centrations, but also to high temperatures and aridity, which enhance water losses through transpiration. Here, the 
environmental factors controlling stomatal behaviour of leaf-level carbon and water exchange were examined across 
the evolutionary continuum from C 3 to C 4 photosynthesis at current (400 umol mol 1 ) and low (280 umol mol 1 ) atmos- 
pheric C0 2 conditions. To this aim, a stomatal optimization model was further developed to describe the evolutionary 
continuum from C 3 to C 4 species within a unified framework. Data on C 3 , three categories of C 3 -C 4 intermediates, 
and C 4 Flaveria species were used to parameterize the stomatal model, including parameters for the marginal water 
use efficiency and the efficiency of the C0 2 -concentrating mechanism (or C 4 pump); these two parameters are inter- 
preted as traits reflecting the stomatal and photosynthetic adjustments during the C 3 to C 4 transformation. Neither 
the marginal water use efficiency nor the C 4 pump strength changed significantly from C 3 to early C 3 -C 4 intermediate 
stages, but both traits significantly increased between early C 3 -C 4 intermediates and the C 4 -like intermediates with 
an operational C 4 cycle. At low C0 2 , net photosynthetic rates showed continuous increases from a C 3 state, across 
the intermediates and towards C 4 photosynthesis, but only C 4 -like intermediates and C 4 species (with an operational 
C 4 cycle) had higher water use efficiencies than C 3 Flaveria. The results demonstrate that both the marginal water use 
efficiency and the C 4 pump strength increase in C 4 Flaveria to improve their photosynthesis and water use efficiency 
compared with C 3 species. These findings emphasize that the advantage of the early intermediate stages is predomi- 
nantly carbon based, not water related. 

Keywords: C 3 -C 4 intermediates, C 3 photosynthesis, leaf gas exchange, photosynthetic model, stomatal conductance, water 
use efficiency. 



Introduction 

While only 3% of the world's terrestrial plant species use the 
C 4 photosynthetic pathway, C 4 species are responsible for 
some 20% of global gross primary productivity (Sage et al, 
2012). The high productivity of C 4 plants is due to their 
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efficient photosynthetic physiology, which includes an addi- 
tional yet spatially separated metabolic cycle, mediated by 
phospho(?«o/pyruvate carboxylase (PEPCase), to the conven- 
tional C 3 Calvin-Benson cycle. This additional cycle results 
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in high C0 2 concentrations around Rubisco, thus suppressing 
the enzyme's oxygenase function and nearly eliminating pho- 
torespiration and the associated carbon and energetic costs. 

C 4 photosynthesis has evolved independently at least 66 
times in lineages throughout the plant kingdom and, in some 
of these lines, there are intermediate species that are neither 
C 3 nor fully C 4 (Sage et al, 2012). Phylogenetic analyses in 
some evolutionary groups where the full range of C 3 , C 3 C 4 
intermediates, and C 4 species are found (such as the genus 
Flaveria) have confirmed that C 3 photosynthesis is the basal 
state and C 4 photosynthesis is more derived; these analyses 
also place photosynthetic intermediate species as evolution- 
ary intermediates to these two photosynthetic types (McKown 
et al, 2005). The C 3 -C 4 intermediates can be classified into 
three categories based on the degree to which they express C 4 
traits: Type I intermediates show refixation of photorespira- 
tory C0 2 by Rubisco in enlarged bundle sheath cells; Type II 
species have increased PEPCase activity and some C 4 cycle 
function; and C 4 -like species have an operational C 4 cycle, but 
have some residual Rubisco activity in their mesophyll cells 
(Edwards and Ku, 1987). 

Despite extensive research, the role of environmental fac- 
tors in driving the evolution of C 4 photosynthetic traits con- 
tinues to draw significant attention (e.g. Osborne and Sack, 
2012; Griffiths et al, 2013). Much work has focused on the 
importance of a drop in atmospheric C0 2 concentrations (c a ) 
from near 1000 ppm to ~400 ppm ~30 million years ago (mya), 
where the predominant benefit from a C0 2 -concentrating 
mechanism would have been enhanced net C0 2 fixation rates 
through suppression of photorespiration (Ehleringer et al, 
1997; Sage, 2004; Christin et al, 2011; Gowik et a/., 2011; 
Sage et al, 2012). While low c a increases photorespiration, 
the effect is even greater when combined with high tempera- 
tures: very low c a conditions, such as those of the last gla- 
cial period (~180 umol mok 1 ) (Luthi et al, 2008), may have 
selected for traits in some C3 species to favour the capture 
and reassimilation of respired and photorespired C0 2 to off- 
set this stress (Busch et al, 2013), but the detrimental effects 
of low c a conditions on C 3 species are further exacerbated 
at warmer temperatures (Campbell et al, 2005). The warm 
regions where C 4 species evolved therefore probably stimu- 
lated photorespiration considerably, but they also drove a 
concomitant increase in transpiration demand (Taylor et al. , 
2012). It has been known for decades that C 4 plants are more 
water use efficient than C 3 species under the same conditions 
(e.g. Rawson et al, 1977; Morison and Gilford, 1983); a spate 
of recent work has highlighted the role of other environmen- 
tal variables that, along with low c a , may have contributed to 
the rise of C 4 photosynthesis, such as dry or saline conditions. 
These recent studies have emphasized the role of C 4 photo- 
synthesis in improving plant water status and preventing 
hydraulic failure in these environments (Osborne and Sack, 
2012; Griffiths et al, 2013). 

Since C 4 species can maintain high photosynthetic rates 
even when stomatal conductance is low compared with their 
C 3 counterparts, it follows that C 4 photosynthesis promotes 
higher water use efficiencies (WUEs) than are found in C 3 
species (e.g. Rawson et al, 1977; Morison and Gifford, 1983; 



Monson, 1989; Huxman and Monson, 2003; Kocacinar et al, 
2008; Osborne and Sack, 2012). This pattern is apparent in 
both forms of WUE: instantaneous WUE (WUE;, defined as 
the leaf carbon gain from net photosynthesis, A Bet , per unit 

d A 

water lost via transpiration, E) and marginal WUE [ A= — — ; 

d E 

see Lloyd and Farquhar (1994); Vogan and Sage (2011); 

Manzoni et al (201 1); note that this definition of X is consist- 
ent with that of Hari et al (1986), but the inverse of the same 
symbol used by Cowan and Farquhar (1977)]. The marginal 
WUE X can also be interpreted as the cost of losing water in 
carbon units. Thus, the higher X of C 4 species implies that 
water loss is more costly for the carbon balance with respect 
to C 3 species, so that C 4 species operate at a relatively low E, 
but at comparable or higher A 1Kt . This finding is consistent 
with C 4 photosynthesis preventing hydraulic failure by means 
of a tight stomatal regulation of water loss (Osborne and 
Sack, 2012). Yet what controls the stomatal behaviour and 
WUE across the evolutionary continuum of C 3 to C 4 species 
remains a subject of research (Vogan and Sage, 2011; Way, 
2012) and frames the scope of this work. 

Recent experiments have shown that instead of a gradual 
improvement in WUE from C 3 species, across the interme- 
diate range and to a full C 4 pathway, the increase in WUE 
resembles a threshold effect: Type I and II intermediates have 
WUEs on a par with C 3 species, while C 4 -like species have a 
high WUE akin to C 4 plants (Kocacinar et al, 2008; Vogan 
and Sage, 2011). The development of the C0 2 -concentrating 
mechanism, which effectively pumps C0 2 from the substoma- 
tal cavity into the chloroplasts where the Calvin-Benson cycle 
occurs, is thought to be the primary mechanism by which 
C 4 plants enhance their WUE ; . The present study therefore 
sought to investigate the connection between the evolutionary 
continuum of C 3 to C 4 photosynthesis and stomatal behav- 
iour (which is a key factor in controlling WUE), by exploring 
the following questions. 

(i) To what degree can stomatal optimization theories describe 
the WUEj patterns in species that have photosynthetic 
characteristics intermediate between C 3 and C 4 species? 

(ii) What are the relationships among the C0 2 -concentrating 
mechanism, X, and WUE ; in these C 3 -C 4 intermediates? 

To address these questions, the genus Flaveria was used as 
a case study, since it contains species with C 3 photosyn- 
thesis, all three intermediate photosynthetic types, and C 4 
photosynthesis (a system previously used by Huxman and 
Monson, 2003; Sage, 2004; McKown and Dengler, 2007; 
Kocacinar et al, 2008; Gowik and Westhoff, 2011; Gowik 
et al, 2011; Vogan and Sage, 2011; and others). Data on 
Flaveria are used to parameterize a stomatal optimization 
model and examine stomatal behaviour across the evolu- 
tionary range of C 3 to C 4 photosynthesis. By using a phy- 
logenetically constrained system, the patterns of changes 
in the model parameters across the C 3 -C 4 photosynthetic 
continuum can be simultaneously explored while minimiz- 
ing evolutionary differences between groups that might oth- 
erwise confound the analysis. 
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Theory 

In the Farquhar et al. (1980) photosynthesis model, A net is 
determined by the minimum of two limitations: the Rubisco 
carboxylation rate (A c ) and ribulose-l,5-bisphosphate 
(RuBP) regeneration rate A h and is commonly expressed as 



A aet = mm(A c ,A J )-R i , 



(1) 



where R d is the daytime respiration rate (see Table 1 for 
symbols and definitions). Rubisco limitation occurs under 
saturating light or at low C0 2 concentrations at the site of 
Rubisco, while RuBP regeneration tends to limit photosyn- 
thesis when c a is high and light levels are low, resulting in a 
limited electron transport rate. The Rubisco-limited assimila- 
tion rate, A c , can be expressed as 



c, + K,, 



(2) 



where F c max is the maximum Rubisco carboxylation rate, c c is 
the C0 2 concentration at the photosynthetic site, T* is the C0 2 
compensation point in the absence of mitochondrial respiration, 
and K cair = K c (1 + O I K 0 ) , with K c and K 0 being the Michaelis- 
Menten constants for Rubisco C0 2 fixation and oxygen inhi- 
bition, and O is the oxygen concentration in the air (21%). 
Conversely, the RuBP-limited assimilation rate is constrained by 
the rate of electron transport, /, and can be expressed as 



/ c-T 



4 c+2T 



(3) 



where the electron transport rate is given by J-a p B m Q, and 
Q is the irradiance, a p is the leaf absorptivity, and e m is the 
maximum photochemical efficiency (Genty et al, 1989). 

To avoid discontinuities in A Bet due to an abrupt transi- 
tion from one limitation to another, the minimum function in 
Equation 1 has often been replaced by a quadratic function, 
at the cost of introducing an additional curvature parameter. 
An alternative approach is to approximate Equation 1 by a 
hyperbolic function (Vico et al, 2013), 



-r* 



(4) 



J 



where K -— andfc 2 -JK cajr /(4K ( , mav ) Such a represen- 
tation ensures that at low c c lk 2 ,A CJ ~A c and at large 
cJ^Acj-Aj. 

When c c I k 2 is approximately unity, both 
Rubisco and RuBP regeneration rates exert comparable limi- 
tations on photosynthesis. Hereafter, this regime is referred to 
as the co-limitation regime. Under C0 2 -limited (or light-satu- 
rated) conditions in which k x - V c max and k 2 -K ah , the optimal 
solution is identical to the one obtained by Katul et al. (2010) 
for non-linear photosynthetic kinetics without light limita- 
tion. Based on Equation 4, net photosynthesis is obtained as 

A nel = Aqj-R^. 

Although there are numerous physiological and ana- 
tomical traits that underlie the development of the C0 2 - 
concentrating mechanism in C 4 plants (e.g. McKown and 



Dengler, 2007; Sage et al, 2012), for modelling purposes, the 
simplest description of the effect of such a pump is to assume 
that the C0 2 concentration at the site where photosynthesis 
occurs is c=r\c w where r| represent the strength of the C0 2 - 
concentrating pump. The value of n encompasses not only 
the development of C 4 biochemistry across the evolutionary 
gradient of species, but also biochemical and anatomical fea- 
tures that affect mesophyll conductance. In C 4 species, r|>l 
(Manzoni et al, 2011); while it is slightly smaller than unity 
in C 3 species due to the need to diffuse C0 2 through the mes- 
ophyll, the lack of specific data on mesophyll conductance 
meant this had to be neglected and r|=l was set for C 3 species. 

In the following, the pump strength, n, is estimated from 
the slope of the A wi {c) curve. Employing this simple descrip- 
tion of the C0 2 -concentrating mechanism results in a simpler 
photosynthesis model than by considering explicitly PEPCase 
kinetics (Collatz et al, 1992; Laisk and Edwards, 2000; von 
Caemmerer, 2000, 2013; Vico and Porporato, 2008), thereby 
allowing data sets collected across different experiments and 
conditions to be compared. Nevertheless, the parameter r| 
can be linked to the kinetics of PEPCase. The C0 2 concentra- 
tion in the stomatal cavity (q) is assumed to be transported by 
PEPCase activity and the shuttling of C 4 acids to the bundle 
sheath (the site of photosynthesis), where the C0 2 concentra- 
tion reaches c bs . When PEPCase kinetics are assumed to be 
linear for illustration (but see von Caemmerer, 2000 for more 
detailed and non-linear models), then 



F P = aq , 



(5) 



where a is the kinetic constant of the process. Setting V p =A net 
(from Equation 4) to guarantee continuity in the C fluxes 
from the stomatal cavity to the site of photosynthesis pro- 
vides an equation to be solved for r\, leading to 



kiT * + k 2 (act + Rj ) 
q (k Y -ac t -R d ) 



(6) 



Equation 6 shows that the pump efficiency r| in principle 
depends on the photosynthetic parameters as well as q. 
However, neglecting R d and assuming r*«q and ac i «A: 1 , 
it can be shown that ?7 = I K > which is a constant at a 
given temperature and light level. Therefore, when respi- 
ration is small and photosynthetic capacity is large, a con- 
stant efficiency r| captures the main effect of the PEPCase 
on the photosynthetic rate. Outside these simplifications, the 
assumption of a constant r| can only be regarded as a first- 
order approximation. 

The combination of the hyperbolic function in Equation 4 
with the simplified description of the C0 2 pumping mechanism 
based on r| (i.e. qpnq) provides a tool to describe C0 2 demand 
within a common framework valid across the C 3 to C 4 evo- 
lutionary continuum. Despite the inherent simplifications, this 
model is in good agreement with earlier, more complex photo- 
synthesis models for C 3 -C 4 intermediates and C 4 species (von 
Caemmerer, 1989; Collatz et al, 1992) (data not shown; for an 
example of model comparison for C 3 species, see fig. 1 in Vico 
et al, 2013), thus lending support to the present approach. 
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Table 1 . Symbols and their definitions used in the paper and model 



Svmboi 


D&finition 


Units 


a 


Ratio of the molecular diffusivities of C0 2 to water vapour 


— 




Fiubisco-limited C0 2 assimilation rate 


Iimol rrr 2 s~ 1 


Aj 


RuBP regeneration-limited net C0 2 assimilation rate 


limol rrr 2 s~ 1 


Aiel 


Net C0 2 assimilation rate 


ixmol rrr 2 s~ 1 


Ca 


Atmospheric C0 2 concentration 


|imol moL 1 




Bundle sheath C0 2 concentration 


limol moL 1 


Co 


Chloroplastic C0 2 concentration 


limol moL 1 


CE 


Carboxylation efficiency 


mol rrr 2 s~ 1 


c. 


Intercellular C0 2 concentration 


(imol mor 1 


D 


Vapour pressure deficit 


kPa 


£ 


Transpiration rate 


mmol rrr 2 s~ 1 




Stomatal conductance 


mmol rrr 2 s~ 1 


J 


Electron transport rate 


limol rrr 2 s~ 1 




Maximum photosynthetic rate of the hyperbolic model (Equation 4) 


limol rrr 2 s~ 1 


/c 2 


Half-saturation constant for the hyperbolic model (Equation 4) 


- 


K= 


Michaelis-Menten constant for Rubisco carboxylation 


limol mor 1 


^cair 


Half-saturated constant of the Rubisco-limited photosynthesis 


limol moL 1 


^cat 


Catalytic constant for Rubisco 


mol mor 1 s~ 1 




Michaelis-Menten constant for Rubisco oxygenation 


mmol moL 11 


o 


Atmospheric oxygen concentration 


mmol mor 1 


Q 


Irradiance 


limol rrr 2 s~ 1 


q 


Temperature coefficient 


- 




Day respiration 


limol rrr 2 s~ 1 




Leaf temperature 


°C 




Maximum carboxylation rate of Rubisco 


limol rrr 2 s~ 1 


^c,max25 


Maximum carboxylation rate of Rubisco at 25 °C 


limol rrr 2 s~ 1 




PEPCase rate 


limol m -2 s~ 1 


WUE 


Water use efficiency 


mmol mor 1 


WUE, 


Instantaneous water use efficiency 


mmol mor 1 


C( 


Kinetic constant (Equations 5 and 6) 


— 


Op 


Leaf absorptivity 


— 


«1.2 


Parameter groups (Equations 8-11) 


— 


P1,2,3 


Parameter groups (Equations 8-11) 




Em 


Maximum photochemical efficiency 




Y 


Parameter groups (Equations 8-11) 




r 


C0 2 compensation point 


limol mor 1 


r* 


C0 2 compensation point in the absence of mitochondrial respiration 


limol mor 1 


n 


C4 pump strength 






Marginal water use efficiency 


mmol mor 1 



The biochemical demand for C0 2 described by A CJ is met 
by C0 2 supplied by the atmosphere via Fickian diffusion at 
a rate given by 

Aet = g s (c a -Ci), (7) 

where g s is the stomatal conductance and c a is the atmos- 
pheric C0 2 concentration. For a given c a , set of environmen- 
tal conditions (such as Q and temperature), and physiological 
properties determining k { and k 2 , the atmospheric supply 
and biochemical demand for A C J constitute two equations 
with three unknowns: g s , c a , and A net . Hence, one addi- 
tional equation is needed to close this system of equations 
mathematically. 

This additional equation can take on the form of an opti- 
mality rule, whereby stomata are assumed to operate so as to 



maximize their carbon gain at a given water loss cost (Cowan 
and Farquhar, 1977; Hari et al, 1986). This hypothesis is 
equivalent to maximizing a Hamiltonian function H-A net - 
XE, where E-ag s D is the leaf transpiration rate (assuming a 
perfectly coupled canopy), and a= 1 .6 is the ratio of the molec- 
ular diffusivities of C0 2 to water vapour. Combining the bio- 
chemical demand with atmospheric supply so as to eliminate 
q, and thereby expressing A C J as a function of g s , inserting 
the outcome into the Hamiltonian, and setting dH ldg s -0, 
leads to a quadratic equation in g s (Vico et al, 2013). Solving 
this equation for g s results in a solution for optimal g s as a 
function of biochemical parameters (n, F c max , AT cair , i? d , and 
T*), environmental conditions (c a and D), and the optimiza- 
tion parameter X. The explicit functional form for optimal 
stomatal conductance is determined from the solution to the 
optimality problem as 
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dH d(A net -AE) optimaUtyconditioB + ^ 

°gs dg s P 3 



where 



^i=-YV(k 2 +a 2 )[k 1 (k 2 -'nc a + 2T*) + a l R d j, (9) 

$ 2 =-ynk 1 (k 2 +T*)(k 2 + a 2 ) 

[k 2 + v(c a -2y)] 2 [k l (T* -^c^+a^], (10) 



p i = yaf(k 2 +a 2 ) 



(11) 



In Equations 9, 10, and 11, a^k^^c^, a 2 =Ti(c a -Y), and y-aXD. 
Therefore, the optimal stomatal conductance depends on X, 
which by using the optimization condition dH ldg s -0 can be 
shown to be equal to the definition of the marginal WUE, i.e. 



dA ni ,f dA„, 



dE dg s 



dE_ 



(12) 



Equation 12 provides a physical interpretation for X, but does 
not give additional information (the optimization condition has 
already been used in Equation 8). Hence, X needs to be determined 
to close the optimization problem mathematically. Although X 
changes as a function of time when soil moisture declines during 
a dry period (Manzoni et al, 2013), under well- watered condi- 
tions or stable moisture levels, X can be considered time-invariant. 

Before applying the proposed model, it is important to 
summarize its key assumptions and simplifications: 

(i) Photosynthetic kinetics are described by a hyperbolic 
function of c { bridging a C0 2 -limited regime (where A mt 
scales linearly with c { ) and a light-limited regime (where 
A nst depends solely on a light level). 

(ii) PEPCase kinetics are described by a single efficiency 
parameter r\, which approximates more complex models 
well (Collatz et al, 1992; von Caemmerer, 2000) when 
respiration terms are small. Also mesophyll resistance is 
neglected, due to a lack of data across these species; this 
assumption implies that the estimated r| could be inflated 
under dry conditions (though these are not the condi- 
tions considered when inferring the marginal WUE). 

(iii) Stomatal conductance is obtained from an optimization 
argument assuming that the marginal WUE is constant — 
a reasonable approximation for experiments under con- 
trolled conditions and stable moisture levels (Manzoni 
et al, 2013). Thus, X is used as a fitting parameter affect- 
ing the stomatal conductance in Equation 8. 

Clearly, these assumptions could be relaxed, thereby improv- 
ing realism. However, relaxing these assumptions reduces the 
ease of interpretation of the derived equations and the abil- 
ity to compare across a wide range of data sets due to more 
required parameters. Once the optimal g s is determined, A c ,i, 
E, and c ; can then be computed. This model allows the quan- 
tification of A net and g s for C 3 , C 4 , and C 3 C 4 intermediate 



species within a common framework and as a function of 
both environmental conditions (air temperature, Q, D, and 
c a ) and species-specific parameters (t|, X, V cm!lx , K C!iiT , and 
T*)). As such, after showing that the modelled response of 
A wX to changes in g s is well captured assuming optimal sto- 
matal behaviour, the model is used to investigate how A net and 
WUEj are altered by changes in n, X, and c a , thus following 
the steps of the hypothesized evolution of C 3 -C 4 intermedi- 
ates and C 4 species from C 3 plants. 



Data availability and model 
parameterization 

To examine the consequences of intermediacy on stomatal 
behaviour, WUE ; , and X, Flaveria species were used where 
gas exchange measurements for C 3 , C 4 , and intermediate spe- 
cies have been previously characterized. The data included C 3 
species (F. cronquistii and F. pringleii), Type I intermediates 
(F. angustifolia, F. chloraefolia, F. pubescens and F. sonoren- 
sis), Type II intermediates (F. floridana and F. ramosissima), 
C 4 -like intermediates (F. brownii, F. palmeri and F. vaginata), 
and C 4 species (F. australasica, F. bidentis, F. kochiana and 
F. trinervia) (see Supplementary Table SI available at JXB 
online). The most recent photosynthetic classification of the 
species (Vogan and Sage, 2011) was employed. Responses of 
g s and A mt to variation in c a , and responses of A net to changes 
in g s were either taken from the literature or digitized from 
published graphs (Monson, 1989; Vogan and Sage, 2011). 
Environmental conditions (Q, D, c a , and leaf temperature) 
in model runs were matched to the measurement conditions 
described for the experimental data. The range in X necessary 
to capture measured responses in gas exchange was explored 
in Flaveria species from all photosynthetic types. 

To parameterize the above model for Flaveria, F c max val- 
ues were derived for Rubisco from 15 Flaveria species that 
spanned C 3 to C 4 photosynthetic types using in vitro meas- 
urements of catalytic constants (or turnover numbers, k cM ) 
(Kubien et al, 2008) and Rubisco site concentrations from 
the same experiment (D. Kubien, personal communica- 
tion) (Table 2). The Michaelis-Menten constants K c and K 0 
for Rubisco were also taken for each Flaveria species from 
Kubien et al, (2008). Rubisco kinetics were adjusted to 30 °C 
to match conditions in the carboxylation efficiency studies 
using correction equations and coefficients from Campbell 
and Norman (1998) (Table 2). K c and K Q were temperature 
adjusted by multiplying their values at 25 °C by exp[^(7' 1 -25)], 
where q is the temperature coefficient for that parameter 
(0.074 for K c and 0.015 for K Q ) and T\ is leaf temperature. 

A . tA T/ F c „„„ A . 25 exp[0.088(r / -25)] 

Pc,max was adjusted as V c , max = — m , n/ „ ttt, — , 

l + exp[0.29(7] -41)] 

where F c max25 is the maximum carboxylation rate at 25 °C. 
The T* values for each of the five photosynthetic types were 
approximated using averaged values of T (the C0 2 compensa- 
tion point) from Flaveria species in Ku et al. (1991), assuming 
that day respiration of mitochondria is small (7? d =0.0 1 5 V c >max ) 
and can be ignored (Table 2). Because light conditions var- 
ied across experiments, the estimation of /(a p e m 0 needed 
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Table 2. Parameter values (based on mean values from experimental data corrected to 30 °C) used for modelling photosynthesis for 
each photosynthetic type of Flaveria 





CE 

(mol nr 2 s~ 1 ) 


r 

(nol mol -1 ) 


1/ 

" c.max 

(nmol m~ 2 s~ 1 ) 


(nmol mol 1 ) 


K 0 

(mmol 
mol 1 ) 


Ad 

(limol m 2 

s 1 ) 


X 

(mmol mol" 1 ) 


1 

(unitless) 


References for 


Sudderth ef al. 


Ku etal. (1991) 


Kubien ef al. 


Kubien ef al. 


Kubien ef al. 


Estimated as 






values 


(2007), citing 




(2008) (and pers. 


(2008) 


(2008) 


0.015 l/ c , max 








Dai etal. (1996) 




comm.) 












C 3 species 


0.11 


61.36 


53.20 


494.9 


575.8 


0.8 


0.826 


1.83 


Type 1 


0.079 


25.45 


73.40 


516.6 


665.7 


1.10 


0.739 


1.73 


Type II 


0.13 


9.20 


66.72 


547.1 


614.6 


1.00 


0.754 


1.93 


C 4 -like 


0.27 


4.93 


34.10 


690.3 


422.3 


0.51 


3.267 


8.83 


C 4 species 


0.47 


3.32 


39.52 


898.7 


1631.8 


0.59 


3.410 


17.60 



CE, carboxylation efficiency; r, C0 2 compensation point; V cmax , In vitro maximum carboxylation rate of Rubisco; K c and K 0 (at 30 °C), Rubisco 
Michaelis-Menten constants for C0 2 and 0 2 , respectively; R a , day respiration rates; X, marginal water use efficiency; tj, C 4 pump strength. 

Data are taken from literature sources as outlined in the text (see also species-specific data points in Fig. 2 and Supplementary Table S1 at 
JXB online); X and r| are calculated values. 



in determining a x = J mda 2 = Jk 2 1 (4V cmax ) requires an 
estimate of the product a p £ m (not their individual values). 
The value for a p was set at 0.8 (based on values for C 3 and 
C 4 species in von Caemmerer, 2000 and Collatz et al, 1992, 
respectively), e m was 0.1 mol mol" 1 (similar to Norman and 
Campbell, 1998; Cheng et al, 2001; Taiz and Zeiger, 2010), 
and Q was set for the irradiance used in individual papers 
being modelled. Hence, a p e m =0.08, resulting in a^O.082 
when RuBP regeneration limits photosynthesis. This esti- 
mate is consistent with conventional values for C 3 species 
(Campbell and Norman, 1998; see table 14.1) though a p e m 
may be more uncertain for C 4 species and C 3 -C 4 intermedi- 
ates. For C3-C4 intermediates, Monson (1989) and Monson 
and Jaeger (1991) report mid-day photosynthetic rates for 
several species including F. floridana between 1 5 umol nT 2 s" 1 
and 45 umol nT 2 s" 1 at light levels ranging from Q- 1 500 umol 
nT 2 s" 1 to 2000 umol nT 2 s -1 . Because Aj ~ J 14 (assuming 
that T* is negligible in Equation 3), it follows that the meas- 
ured A are consistent with the estimate J-0.08Q, which gives 
an A net of 40 umol nT 2 s" 1 . This evidence supports the assump- 
tion that a p £ m is stable across photosynthetic types. A more 
rigorous parameterization would require direct observations 
of a p £ m or reliable V cm:ix -J m ^ relationships across the C 3 C 4 
continuum. 

Carboxylation efficiencies [CEs; i.e, the initial slope of the 
v4 net (ci) curve] measured under saturating light for Flaveria 
species were taken from Krall et al. (1991) and Sudderth et al. 
(2007), citing Dai et al. (1996). Based on Equation 2, at low q 
the slope of the ^4 net (Ci) is approximately CE = V c max 7] I K cair , 
so that ?7 = CEK mir /V c max . Therefore, knowledge of T, K ail , 
K c , and K a allowed an estimate of n for each Flaveria species. 

Finally, the Vogan and Sage (2011) gas exchange data set 
was used to infer how X changes across the evolutionary path- 
way from C 3 to C 4 species. In that study, a range of g s and 
A net values was obtained by altering the nitrogen availability 
for individuals of all photosynthetic types considered, while 
water was amply supplied, so X can be considered time-invar- 
iant. As a consequence of different nutrient availability, a 
range of photosynthetic capacities and respiration rates were 



obtained. Since there is no way of knowing these biochemi- 
cal parameters, a simplified but more robust approach to 
estimate X was adopted that only requires gas exchange rates 
and photosynthetic type-averaged T and r| (assuming X is 
substantially unaltered by nutrient availability). For this step, 
instead of using the definition (Equation 12), which requires 
knowledge of all the photosynthetic parameters, the stoma- 
tal optimization model was simplified by selecting light-satu- 
rated conditions, so that R d ~0 and the photosynthesis model 
is approximately linear. Following these simplifications, it can 
be shown that A net = g s JaDA,(c a —T* It]) (Manzoni et al. , 
201 1), which allows estimating X through a linear least square 
regression of A mt versus g s constrained through the origin for 
each photosynthetic type. In previous works on different spe- 
cies, this approach to estimate X was compared with results 
obtained without these simplifications. Such comparison 
showed that the differences between the two approaches was 
rarely more than 20% across a wide range of environmental 
and physiological conditions (see fig. 4 in Katul et al, 2010), 
which is in the range of experimental variability [e.g. a mean 
standard deviation of 16% in light-saturated A net estimates 
across individuals in a range of C 3 , C 3 C 4 intermediates, and 
C 4 species (Vogan et al, 2007)]. 

Gas exchange rates were also simulated under altered 
atmospheric C0 2 concentrations. In this analysis, all bio- 
chemical parameters were maintained constant, but the pos- 
sibility was considered that the marginal WUE increases 
linearly with C0 2 concentrations (Manzoni et al, 2011). 
Simulations with constant X estimated as described above 
were thus compared with simulations with X(c a )=^ 40 o(c a /400). 
Including C0 2 effects allows the robustness of the results to 
changes in X to be tested. 

Results and Discussion 

Recent work has stimulated new interest in the role that tran- 
spiration demands may have played in the evolution of C 4 
photosynthesis and traits associated with the C 4 syndrome 
(Taylor et al, 2011, 2012; Osborne and Sack, 2012; Griffiths 
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et al, 2013). These studies have emphasized that C 4 photo- 
synthesis not only benefits the carbon economy of a plant, 
but also has important implications for hydraulic traits, 
drought tolerance, and water use patterns, benefits that are 
maintained or enhanced when C 4 plants are exposed to the 
low c fl conditions where C 4 photosynthesis evolved (Ripley 
et al. , 2013). Here, X is used as an 'index' of the cost of losing 
water in terms of carbon, and its variation along the evolu- 
tionary gradient from C 3 to C 4 photosynthesis is investigated. 

Combining a stomatal optimization approach with meas- 
ured biochemical parameters, realistic mean ^4 net (Cj) curves 
for each photosynthetic pathway in Flaveria were computed 
(Fig. 1A; Table 2). The corresponding estimated n values 
are reported in Fig. IB. In the optimality model, recall that 
the parameter n represents an overall pump strength for the 
carbon-concentrating mechanism, which might naively have 
been expected to increase gradually from C 3 towards C 4 
plants. Instead, the analysis here suggests that n was relatively 
stable and similar to that for C 3 species (rpl or slightly above 
1 due to unavoidable errors in the estimation) until reaching 
C 4 -like species. The relatively constant r| between C 3 , Type I, 
and Type II Flaveria species occurred despite there being an 
increase in the initial slope of the ^4 net (cj) curve (e.g. the car- 
boxylation efficiency) across these groups. Instead of being 
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Fig. 1 . (A) Responses of net C0 2 assimilation rates (A net ) to increases in 
intercellular C0 2 concentration (q), relativized to maximum A net (=/4 max ) for 
each photosynthetic type, as commonly presented in the literature. (B) 
Estimated r| for each photosynthetic type. Means ±SE, n indicated at the 
top, dotted line indicates r|=1 . C 3 species, purple circles and solid line; 
Type I species, blue diamonds and dashed line; Type II species, green 
triangles and dotted line; C 4 -like species, yellow inverted triangles and 
dashed-dotted line; C 4 species, red squares and solid line. 



attributed to r|, the steeper initial slopes in the Type I and 
Type II intermediates in comparison with the C 3 species were 
caused by higher F c max values for Rubisco based on in vitro 
assays of the enzyme kinetic parameters (Table 2; Fig. 2A), 
consistent with positive selection on Rubisco across the C 3 
to C 4 gradient (Kapralov et al., 2011). Thus, there was no 
increase in n until the C 4 -like species were reached; at this 
point, T] values were about half-way between those of the 
full C 3 and C 4 photosynthetic groups. The greater variation 
in t] estimates in species closer to the C 4 end of the spec- 
trum is therefore probably due to the greater range of pump 
strengths possible as the carbon-concentrating mechanism is 
established and to the species-level diversity in V Q max values 
(Supplementary Table SI at JXB online). While there were 
sharp changes in in vitro Rubisco V^ m!lx between C 3 species 
and the Type I and Type II intermediates, the change in the 
of Rubisco across the photosynthetic groups was more 
gradual until reaching the C 4 -like species (Fig. 2B), implying 
that these enzyme kinetic traits are not necessarily linked. The 
T* dropped sharply as r| increased slightly above a value of 1 
and then flattened (Fig. 2C). 

While it might, a priori, seem reasonable to expect a 
steady increase in WUE; from C 3 species through the inter- 
mediate Flaveria species and to full C 4 plants, this was not 
borne out by the data, in agreement with published findings. 
Huxman and Monson (2003) showed that the WUEjS of 
C 3 -C 4 intermediates were similar to C 3 WUEj values, while 
C 4 WUEj values were considerably higher. Vogan and Sage 
(2011) also found no evidence for a gradual transition in the 
slope between A net and g s in C 3 C 4 intermediates, but rather 
a sharp increase between Type II intermediates and C 4 -like 
intermediates. [Note that WUE; is a proxy for X if a linear 
^net^D curve is assumed.] In our re-analysis of the Vogan and 
Sage (201 1) data set, the reported gas exchange data could be 
readily described with the optimality model for C 3 , and Type 
I and Type II intermediates using the estimated changes in 
n (as in Fig. IB), but without significant changes in X across 
photosynthetic pathways (Fig. 3). This result implies that 
there is little change in the relationship between carbon and 
water from that of a C 3 species in these early intermediate 
steps. However, in C 4 -like intermediates, a strong C 4 pump 
(i.e. r|=8) was accompanied by a quadrupling of X compared 
with that used to characterize the data for Type I and II inter- 
mediates. Thus, the stomatal optimization approach could be 
successfully used to capture key changes in the measured rela- 
tionship between A net and g s across the C 3 -C 4 spectrum using 
the estimated r\ values, but only when the marginal WUE of 
the C 4 -like and C 4 species was modelled to be 4-fold greater 
than that of the C 3 species (Fig. 3). This increase in X unam- 
biguously indicates a higher carbon cost for losing water in 
the C 4 -like and C 4 species. In the optimization model, the 
long-term c,/c a dictates X, because Xoo( l-c/c a ) 2 . Therefore, 
the increase in X is generated by a decline in c, (where c a is 
assumed to be 400 umol mok 1 ). This finding suggests that 
the increase in C 4 WUE ; values, the increase in X, and the 
decline in c ; (as expected by the presence of a C 4 pump) are 
all interconnected and predicted from the proposed stomatal 
optimization model. 
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Fig. 2. Relationships between the pump strength of the C0 2 - 
concentrating mechanism (t|) and: (A) maximum carboxylation capacity 
of Rubisco (\/ CjrTlax ); (B) K cajr = K c (1 + OIK 0 ) ; and (C) the C0 2 
compensation point in the absence of mitochondrial respiration (r*) for 
C 3 , C 4 , and C 3 -C 4 intermediate species. Solid lines are fit to data; vertical 
dotted lines indicate n=1 . C 3 species, purple circles; Type I species, blue 
diamonds; Type II species, green triangles; C 4 -like species, yellow inverted 
triangles; C 4 species, red squares. 



Since C 4 photosynthesis evolved under low c a , with the 
transition in Flaveria occurring within the last 3 million years 
(Sage et al, 2012), the effect of varying n and X on A net and 



WUEj was further investigated under both current (400 umol 
mok 1 ) and low C0 2 concentrations (280 umol mol -1 ; Figs 
4, 5), allowing F cmax , K ail , and T* to vary along with n as 
per the relationships in Figs 1 and 2, but keeping X. constant. 
In the results from both current and low c a levels, increases 
in r| initially induce a sharp increase in A net , as more C0 2 is 
concentrated around Rubisco, with a diminishing response 
above a certain r| (rplO for 400 umol moh 1 C0 2 ; Fig. 4A). 
Moreover, higher values of X decrease A net in both environ- 
ments, so that the slight changes in r) and X in Type I and II 
species compared with C 3 plants generate a similar A ne( in the 
three groups (Fig. 4A). Increased c a (from 280 umol m~ 2 s _I 
to 400 umol moF 1 ) elevated A mt in C 3 species by 80% due to 
greater substrate availability (Fig. 4B). Compared with a C 3 
Flaveria at these low c a , C 3 -C 4 intermediates also have higher 
A net in modern C0 2 concentrations, with a gradual increase 
in the stimulation of A nst with respect to C 3 values (Fig. 4B). 
In contrast to the A net results, increases in r\ have little impact 
on WUEj when X is small, namely from C 3 species to Type 
I or II species (shown for 400 umol moF 1 C0 2 in Fig. 5A). 
Moreover, while increases in c a have increased WUEj of 
C 3 species by 30%, there is no gradual rise in WUE ; across 
the gradient of photosynthetic types, as there was with A net 
(Fig. 5B). Instead, compared with a C 3 Flaveria at 280 umol 
moh 1 C0 2 , Type I and Type II intermediates have a similar 
30% stimulation in WUE ; , while C 4 -like and C 4 species show 
a more than tripling of their WUE; stimulation at modern 
C0 2 levels (Fig. 5B). 

If the marginal WUE is assumed to increase with atmos- 
pheric C0 2 (e.g. Katul et al, 2010; Manzoni et al, 201 1), the 
predicted g s at c a =280 umol moF 1 increases. As a consequence, 
photosynthesis also increases and the ratios of net photosyn- 
thesis at c a =400 umol m~ 2 s -1 and 280 umol moF 1 therefore 
decrease (Fig. 4B). Because the positive effect of changes in 
X is larger on transpiration than on net photosynthesis, the 
WUE! at the lower c a decreases. As a result, the ratio of WUEj 
at current and low C0 2 concentrations is higher than when 
assuming a constant X, (Fig. 5B). 

The modelled changes in leaf-level performance between 
photosynthetic groups under low c a are shown in Fig. 6. 
This figure quantifies the advantages of the intermediate and 
C 4 species over the basal C 3 state. At low c a , the estimated 
changes in n and X for intermediate species provide a con- 
tinuous, smooth gradient of increasing carbon gain, over a 
C 3 Flaveria species (Fig. 6A). This trend is robust to changes 
in the X(c a ) relationship, as indicated by minor differences 
between filled and open symbols. A Type I intermediate has 
a 15% higher A net than a C 3 species, which could provide a 
competitive edge to the intermediate in a low C0 2 environ- 
ment; a similar jump in A net is seen for each photosynthetic 
group along the evolutionary trajectory, in agreement with a 
recently proposed smoothly increasing fitness landscape for 
C 4 evolution (Heckmann et al. , 2013). However, the same pat- 
tern is not apparent in the WUEj results (Fig. 6B). There is 
no difference in the WUE ; estimated at 280 umol mok 1 C0 2 
between C 3 , Type I and Type II Flaveria species. Instead, sig- 
nificant increases in WUEj are only achieved in C 4 -like and 
C 4 species, implying that the driving force for the initial steps 
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Fig. 3. Relationships between stomatal conductance to C0 2 (g s ) and net C0 2 assimilation rate (A nBt ) in Flaveria species across the C 3 to C 4 
photosynthetic range. Data points are from Vogan and Sage (201 1); lines are obtained by analytical least-square fitting of the water use efficiency X, 
employing a linearized version of the stomatal optimization model (Manzoni ef al. , 201 1 ) for analytical tractability. 
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Fig. 4. (A) Modelled relationships between net C0 2 assimilation rate (/\ nat ), 
marginal water use efficiency (X), and the C0 2 -concentrating pump strength 
(t|) modelled at current C0 2 concentrations (400 ^mol mol" 1 ); V c max , K cA , 
and r* vary with iq according to the relationships in Fig. 2; vapour pressure 
deficit (D) was set to 1 .5 kPa, leaf temperature to 30 °C, Q to 1500 nmol 
m" 2 s" 1 . Mean values of X and r| for each of the five photosynthetic types 
are indicated on the surface. (B) The ratio of A net at current atmospheric 
C0 2 levels versus A nst of C 3 Flaveria at low atmospheric C0 2 concentrations 
(280 (jmol mol" 1 ) (/4 net , 400 /4 ne , t C3 , 2 so) for each photosynthetic group; means 
±SE across species; filled symbols refer to constant X, open symbols to X 
increasing linearly with c a ; the dashed-dotted line indicates a ratio of 1 . C 3 
species, purple circle; Type I species, blue diamond; Type II species, green 
triangle; C 4 -like species, yellow inverted triangle; C 4 species, red square. 
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Fig. 5. Modelled relationships between instantaneous water use 
efficiency (WUE;, the ratio of A net to £), marginal WUE (X), and the C0 2 - 
concentrating pump strength (r|) modelled at current C0 2 concentrations 
(400 \imo\ mol" 1 ); l/ c max , K cair , and r* vary with r| according to the 
relationships in Fig. 2; vapour pressure deficit (D) was set to 1 .5 kPa, leaf 
temperature to 30 °C, Q to 1500 nmol m" 2 s" 1 . Mean values of X and n 
for each of the five photosynthetic types are indicated on the surface. (B) 
The ratio of WUE, at current atmospheric C0 2 levels versus the WUE, of 
a C 3 Flaveria at low atmospheric 0O 2 concentrations (280 (imol mol" 1 ) 
(WUE i400 /WUE C3 i280 ) for each photosynthetic group; means ±SE across 
species; filled symbols refer to constant X, open symbols to X increasing 
linearly with c a ; the dashed-dotted line indicates a ratio of 1 . C 3 species, 
purple circle; Type I species, blue diamond; Type II species, green 
triangle; C 4 -like species, yellow inverted triangle; C 4 species, red square. 



3692 | Wayefa/. 




Type I 



Type II 



C4-like 




Type I 



Type II 



C4-like 



Fig. 6. Comparison of (A) modelled photosynthetic rates and (B) modelled 
WUEi among Flaveria species with different photosynthetic types at 
c a =280 umol mol" 1 , expressed as ratios over the mean A net and WUE, for 
C 3 species at c a =280 umol mol" 1 . Symbols represent means ±SE across 
species (for fixed C 3 A net and WUEi values); filled symbols refer to constant 
X, open symbols to X increasing linearly with c a ; the dashed-dotted line 
indicates a ratio of 1 . Other parameters are as in Figs 4 and 5. 

towards C 4 photosynthesis in this group was carbon based 
and not related to increasing WUE^ 

Many of the features considered to pre-adapt a group 
to evolve C 4 photosynthesis are related to leaf hydraulics, 
including increased vein density and enlarged bundle sheath 
cell size (McKown et al, 2005; McKown and Dengler, 2007; 
Osborne and Sack, 2012; Sage et al, 2012; Griffiths et al, 
2013). Stomatal anatomy also evolves along the transition 
from C 3 to C 4 photosynthesis, with C 4 species having lower 
maximum stomatal conductance (due to either lower sto- 
matal density or smaller stomatal size) than C 3 congeners 
(Taylor et al, 2012). While changes in whole-plant physiology 
are outside the scope of this work, these findings have stimu- 
lated interest in the role of plant water relations in the evolu- 
tion of C 4 photosynthesis. The results here indicate that while 
there is a gradual increase in carbon gain across the range 
from C 3 to C 4 , there is no corresponding transition in either 
WUEi or "k. Rather, increases in leaf-level WUE ; are only seen 
between Type II intermediacy and C 4 -like species (as noted 
by Kocacinar et al, 2008; Vogan and Sage, 2011). However, 
this transition is accompanied by a rise in X, indicating that a 
higher carbon cost is being incurred for water loss in C 4 -like 
and C 4 species than in C 3 , or Type I or Type II intermediate 



species of Flaveria. This corresponds to the coordinated set 
of changes to the hydraulic architecture of Flaveria spe- 
cies, including lower leaf specific hydraulic conductivity and 
greater cavitation resistance in C 4 and C 4 -like than C 3 species 
(Kocacinar et al , 2008), emphasizing the importance of the 
transition from having a functional C 4 cycle for both the car- 
bon and water economies of the plant. 



Conclusions 

Using a stomatal optimization approach, the full range of C 3 , 
C 3 C 4 intermediates, and C 4 gas exchange could be realistically 
modelled with the addition of a C 4 pump strength parameter 
T], describing the effects of the C 4 carbon-concentrating mech- 
anism. The results here showed that, to capture the patterns 
apparent in measured gas exchange data, the carbon-based 
cost of losing water (k) between C 3 , and Type I and Type II 
intermediates could be maintained constant, but X had to be 
quadrupled to model C 4 -like and C 4 Flaveria (at least within the 
confines of the optimality assumption of stomata). When leaf- 
level fluxes were modelled at low C0 2 , there was no evidence 
for a greater WUE ; in the C 3 -C 4 intermediates (compared with 
a C 3 Flaveria) until they developed a full C 4 cycle. However, 
the model results suggest a steady increase in net carbon fixa- 
tion rates across the C 3 to C 4 photosynthetic range. While this 
implies that carbon, not water, was the main driving pressure 
for the early steps of C 4 evolution in this genus, the increase in X 
indicates that there was a fundamental shift over the evolution 
of C 4 photosynthesis between the relative costs of carbon and 
water, resulting in higher carbon costs of water losses. 



Supplementary data 

Supplementary data are available at JXB online. 
Table SI. Parameter table. 



Acknowledgements 



The authors would like to thank the editor, Rowan Sage, and two anonymous 
reviewers for their constructive comments on the paper. This work was sup- 
ported by grants from the US Department of Agriculture, Agriculture and 
Food Research Initiative (#2011-67003-30222) and the US Department of 
Energy, Terrestrial Ecosystem Sciences (#ll-DE-SC-0006967) to DAW and 
GGK, and from the US-Israeli Bi-national Science Foundation (#2010320) 
and the Natural Sciences and Engineering Research Council of Canada to 
DAW. GV gratefully acknowledges the support of the project AgResource - 
Resource Allocation in Agriculture', from the Faculty of Natural Resources 
and Agricultural Sciences, Swedish University of Agricultural Sciences 
(SLU); SM was supported through an excellence grant from the Faculty of 
Natural Resources and Agricultural Sciences and the vice-chancellor of the 
SLU. 



References 

Busch FA, Sage TL, Cousins AB, Sage RF. 2013. C3 plants enhance 
rates of photosynthesis by reassimilating photorespired and respired C0 2 . 
Plant, Cell and Environment 36, 200-21 2. 

Campbell GS, Norman JM. 1998. An introduction to environmental 
biophysics , 2nd edn. Berlin: Springer. 

Campbell CD, Sage RF, Kocacinar F, Way DA. 2005. Estimation of the 
whole-plant C0 2 compensation point of tobacco (Nicotiana tabacum L). 
Global Change Biology 11,1 956-1 967. 



Stomatal optimization and the transition from C 3 to C 4 photosynthesis | 3693 



Cheng L, Fuchigami LH, Breen PJ. 2001 . The relationship between 
photosystem II efficiency and quantum yield for C0 2 assimilation is not 
affected by nitrogen content in apple leaves. Journal of Experimental 
Botany 52, 1865-1872. 

Christin PA, Osborne CP, Sage RF, Arakaki M, Edwards EJ. 201 1 . 

C 4 eudicots are not younger than C 4 monocots. Journal of Experimental 
Botany 62, 3171-3181. 

Collatz, GJ, Ribas-Carbo M, Berry JA. 1992. Coupled photosynthesis- 
stomatal conductance model for leaves of C4 plants. Australian Journal of 
Plant Physiology 1 9, 51 9-538. 

Cowan IR, Farquhar GD. 1977. Stomatal function in relation to leaf 
metabolism and environment. Symposia of the Society for Experimental 
Biology 31, 471-505. 

Dai Z, Ku MSB, Edwards GE. 1 996. Oxygen sensitivity of 
photosynthesis and photorespiration in different photosynthetic types in 
the genus Flaveria. Planta 198, 563-571 . 

Edwards GE, Ku MSB. 1987. Biochemistry of C3-C4 intermediates. In: 
Strumpf PK, Conn EE, eds.TAie biochemistry of plants. London: Academic 
Press, 275-325. 

Ehleringer JR, Cerling TE, Helliker BR. 1997. C4 photosynthesis, 
atmospheric C0 2 , and climate. Oecologia 112, 285-299. 

Farquhar G, von Caemmerer S, Berry JA. 1 980. A biochemical model of 
photosynthetic C0 2 assimilation in leaves of C3 species. Planta 149, 78-90. 

Genty B, Briantais JM, Baker NR. 1989. The relationship between 
the quantum yield of photosynthetic electron transport and quenching of 
chlorophyll fluorescence. Biochimica et Biophysica Acta 990, 87-92. 

Gowik U, Brautigam A, Weber KL, Weber APM, Westhoff P. 201 1 . 

Evolution of C4 photosynthesis in the genus Flaveria— how many and 
which genes does it take to make C4? The Plant Cell 23, 2087-21 05. 

Gowik U, Westhoff P. 201 1 . The path from C 3 to C 4 photosynthesis. 
Plant Physiology 155, 56-63. 

Griffiths H, Weller G, Toy LF, Dennis RJ. 2013. You're so vein: bundle 
sheath physiology, phylogeny and evolution in C3 and C4 plants. Plant, 
Cell and Environment 36, 246-261 . 

Hari P, Makela A, Korpilahti E, Holmberg M. 1986. Optimal control of 
gas exchange. Tree Physiology 2, 1 69-1 75. 

Heckmann D, Schulze S, Denton A, Gowik U, Westhoff P, Weber 
APM, Lercher MJ. 2013. Predicting C4 photosynthesis evolution: 
modular, individually adaptive steps on a Mount Fuji fitness landscape. Cell 
153, 1579-1588. 

Huxman TE, Monson RK. 2003. Stomatal responses of C3, C3-C4 
and C4 Flaveria species to light and intercellular C0 2 concentration: 
implications for the evolution of stomatal behavior. Plant, Cell and 
Environment 26, 313-322. 

Kapralov MV, Kubien DS, Andersson I, Filatov DA. 201 1 . Changes 
in Rubisco kinetics during the evolution of C4 photosynthesis in Flaveria 
(Asteraceae) are associated with positive selection on genes encoding the 
enzyme. Molecular Biology and Evolution 28, 1 491-1 503. 

Katul G, Manzoni S, Palmroth S, Oren R. 201 0. A stomatal optimization 
theory to describe the effects of atmospheric C0 2 on leaf photosynthesis 
and transpiration. Annals of Botany 105, 431-442. 

Kocacinar F, McKown AD, Sage TL, Sage RF. 2008. Photosynthetic 
pathway influences xylem structure and function in Flaveria (Asteraceae). 
Plant, Cell and Environment 31, 1363-1376. 

Krall JP, Edwards GE, Ku MSB. 1991 . Quantum yield of photosystem II 
and efficiency of C0 2 fixation in Flaveria (Asteraceae) species under varying 
light and C0 2 . Australian Journal of Plant Physiology 18, 369-383. 

Ku MSB, Wu J, Dai Z, Scott RA, Chu C, Edwards GE. 1991 
Photosynthetic and photorespiratory characteristics of Flaveria species. 
Plant Physiology 96, 51 8-528. 

Kubien DS, Whitney SM, Moore PV, Jesson LK. 2008. The biochemistry 
of Rubisco in Flaveria. Journal of Experimental Botany 59, 1 767-1 777. 

Laisk A, Edwards GE. 2000. A mathematical model of C4 
photosynthesis: the mechanism of concentrating C0 2 in NADP-malic 
enzyme type species. Photosynthesis Research 66, 199-224. 

Lloyd J, Farquhar GD. 1994. C-13 discrimination during C0 2 assimilation 
by the terrestrial biosphere. Oecologia 99, 201-215. 

Liithi D, Le Floch M, Bereiter B, ef al. 2008. High-resolution carbon 
dioxide concentration record 650,000-800,000 years before present. 
Nature 453, 379-382. 



Manzoni S, Vico G, Katul G, Fay PA, Polley W, Palmroth S, 
Porporato A. 201 1 . Optimizing stomatal conductance for maximum 
carbon gain under water stress: a meta-analysis across plant functional 
types and climates. Functional Ecology 25, 456-467. 

Manzoni S, Vico G, Palmroth S, Porporato A, Katul G. 201 3. 

Optimization of stomatal conductance for maximum carbon gain under 
dynamic soil moisture. Advances in Water Resources 62, 90-105. 

McKown AD, Dengler NG. 2007. Key innovations in the evolution of 
Kranz anatomy and C 4 vein pattern in Flaveria (Asteraceae). American 
Journal of Botany 94, 382-399. 

McKown AD, Moncalvo JM, Dengler NG. 2005. Phylogeny of Flaveria 
(Asteraceae) and inference of C4 photosynthesis evolution. American 
Journal of Botany 92, 1 91 1 -1 928. 

Monson RK. 1989. The relative contribution of reduced photorespiration, 
and improved water- and nitrogen-use efficiencies, to the advantages of 
C 3 -C 4 intermediate photosynthesis in Flaveria. Oecologia 80, 215-221 . 

Monson RK, Jaeger CH. 1 991 . Photosynthetic characteristics of the 
C3-C4 intermediate Flaveria floridana J.R. Johnson measured in its natural 
habitat. Evidence of the advantages to C3-C4 photosynthesis at high leaf 
temperatures. American Journal of Botany 78, 795-800. 

Morison JIL, Gifford RM. 1983. Stomatal sensitivity to carbon dioxide 
and humidity: a comparison of two C3 and two C4 grasses. Plant 
Physiology 71 , 789-796. 

Osborne CP, Sack L. 2012. Evolution of C4 plants: a new hypothesis for 
an interaction of C0 2 and water relations mediated by plant hydraulics. 
Philosophical Transactions of the Royal Society B: Biological Sciences 
367, 583-600. 

Rawson HM, Begg JE, Woodward RG. 1977. The effect of atmospheric 
humidity on photosynthesis, transpiration and water use efficiency of 
leaves of several plant species. Planta 134, 5-10. 

Ripley BS, Cunniff J, Osborne CP. 2013. Photosynthetic acclimation 
and resource use by the C3 and C4 subspecies of Alloteropsis semialata 
in low C0 2 atmospheres. Global Change Biology 19, 900-910. 

Sage RF. 2004. The evolution of C 4 photosynthesis. New Phytologist 161, 
341-370. 

Sage RF, Sage TL, Kocacinar F. 2012. Photorespiration and the 
evolution of C4 photosynthesis. Annual Review of Plant Biology 63, 
1 9-47. 

Sudderth EA, Muhaidat RM, McKown AD, Kocacinar F, Sage RF. 

2007. Leaf anatomy, gas exchange and photosynthetic enzyme activity in 
Flaveria kochiana. Functional Plant Biology 34, 1 1 8-1 29. 

Taiz L, Zeiger E. 2010. Plant physiology , 5th edn. Sunderland, MA: 
Sinauer Associates. 

Taylor SH, Ripley BS, Woodward Fl, Osborne CP. 201 1 . Drought 
limitation of photosynthesis differs between C3 and C4 grass species in a 
comparative experiment. Plant, Cell and Environment 34, 65-75. 

Taylor SH, Franks P, Hulme S, Spiggs E, Christin P-A, Edwards E, 
Woodward I, Osborne C. 2012. Photosynthetic pathway and ecological 
adaptation explain stomatal trait diversity amongst grasses. New 
Phytologist 193, 387-396. 

Vico G, Manzoni S, Palmroth S, Weih M, Katul G. 2013. A perspective 
on optimal leaf stomatal conductance under C0 2 and light co-limitations. 
Agricultural and Forest Meteorology 1 82-1 83, 1 91 -1 99. 

Vico G, Porporato A. 2008. Modelling C3 and C4 photosynthesis under 
water-stressed conditions. Plant and Soil 313, 187-203. 

Vogan PJ, Frohlich MW, Sage RF. 2007. The functional significance of 
C3-C4 intermediate traits in Heliotropium L. (Boraginaceae): gas exchange 
perspectives. Plant, Cell and Environment 30, 1337-1345. 

Vogan PJ, Sage RF. 201 1 . Water-use efficiency and nitrogen-use 
efficiency of C3-C4 intermediate species of Flaveria Juss. (Asteraceae). 
Plant, Cell and Environment 34, 1415-1430. 

von Caemmerer S. 1 989. A model of photosynthetic C0 2 assimilation 
and carbon-isotope discrimination in leaves of certain C3-C4 
intermediates. Planta 178, 463-474. 

von Caemmerer S. 2000. Biochemical models of leaf photosynthesis. 
Collingwood, Victoria: CSIRO. 

von Caemmerer S. 201 3. Steady state models of photosynthesis. Plant, 
Cell and Environment 36, 1 61 7-1 630. 

Way DA. 2012. What lies between: the evolution of stomatal traits on the 
road to C 4 photosynthesis. New Phytologist 193, 291-293. 



